A1 = readmatrix('10V_alongP.txt')

x1 = A1(:,1)*2*pi/360
y1 = A1(:,2)*1e12*1.2

a =  0.1402*1.2
b =  -0.001853*1.2
d =  1.045*1.2
e =  0.005713*1.2

f1 = a*sin(4*x1)+b*cos(4*x1)+d+e*sin(2*x1)
f2 = e*sin(2*x1)

hold on

%plot(x1,f1,'Linewidth',2,'Color',[0 0 0.8])
plot(x1,f2,'Linewidth',2,'Color',[0.8 0 0])
%scatter(x1,y1,10,[0 0 0.8],'filled')
xlabel('Polarization')
ylabel('Photocurrent [pA]')
xlim([pi, 2*pi])
ylim([-0.15,0.15])
xticks(pi:pi/4:2*pi)
yticks(-0.1:0.1:0.1)
xticklabels({})
yticklabels({})
%set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 200, 80])